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SUMMARY 


A method is suggested for predicting the stability of automati- 
cally controlled aircraft by a comparison of calculated frequency- 
response curves for the aircraft and experimentally determined 
frequency— response curves for the automatic pilot. The method is 
applied only to stabilization in roll. The method is expected to be 
useful as a means of establishing the specifications of the performance 
required of the automatic control device for pilotless aircraft 
designed as missiles. 


INTRODUCTION 


Experience has shown that the provision of automatic stabilization 
for small pilotless aircraft designed as missiles is extremely diffi- 
cult. The difficulty is a result of the high-frequency oscillations of 
small-size aircraft that require rapid control movements and small time 
lags, characteristics which are difficult to obtain, particu lar ly when 
the space available for the control servomotors and intelligence units 
is considered. In an unpublished analysis made at the Langley 
Aeronautical Laboratory of the NACA, the problem of dete r mi n ing the 
stability of an automatically controlled aircraft with lag in the 
control system was analyzed theoretically by ass uming a simplified 
equation for the control motion, this equation being obtained from the 
knowledge of the behavior of the automatic pilot. Because of the 
irregular response characteristics often found in automatic pilots, 
however, the control motion is difficult to represent mathematically 
and, hence, the simplified equations of the control were found to be 
inadequate for the analysis. 



2 


NACA TN 1901 


The present paper suggests a method for predicting the stability 
of an aircraft based on the experimental determination of the character- 
istics of its automatic pilot. The procedure consists essentially in 
calculating the control motion required to maintain a continuous 
sinusoidal motion of unit amplitude for the degree of freedom being 
inspected. The motion of the control is obtained for a range of 
frequencies; the phase angle of the control motion and the ratio of 
amplitude of control motion to airplane motion are plotted as a function 
of frequency. Similar curves are established for the autopilot by 
oscillating it and recording the control motion. The two sets of data 
are then compared to determine whether the airplane will be stable tinder 
control of the automatic pilot. The method is developed in detail only 
for stabilization in roll. It may be used by the airplane designer for 
either determining the suitability of an. existing automatic pilot for a 
particular application or specifying the characteristics of the auto- 
matic pilot needed for the application. 


SYMBOIS 


m 

1 

5 
b 

P 

6 



D 

(0 


mass of airplane, slugs 

radius of gyration of airplane about longitudinal axis, feet 
dynamic pressure, pounds per square foot 
wing area, square feet 
wing span, feet 

rolling-moment coefficient (Bolling mament/qSb) 
angle of bank, radians 

angular velocity in bank, radians per second (d0/dt) 
deflection of aileron, radians 

rate of change of rolling-moment coefficient with a ng u l ar velocity 
in bank, per radian (dCj/dp) 

rate of change of with 8, per radian (dCj/^) 

differential operator (d/dt) 

pingnlR-r frequency, radians per second 
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9 phase angle (positive value means lead of 5 ahead of fi) 
ftma-ir ma ximu m amplitude of 

K control-amplitude ratio (ratio of control deflection to airplane 
displacement) 

l lag in seconds between signal for control and its actual motion 
t time, seconds 

p. real part of root of stability equation 

T ]y 2 for oscillation to damp to one— half its amplitude, seconds 

6(t) control motion as a function of time 
T period of oscillation, seconds 


DETERMINATION OF CONDITIONS FOR NEUTRAL STABILITY 


The method of determining the conditions for neutral stability is 
illustrated in figure 1. The calculated phase angle of the control 
motion and the calculated ratio of the amplitude of control motion to 
airplane motion are plotted against angular frequency as shown by the 
solid— line curves. The upper dashed curve is a plot of the experi- 
mental ratio of the amplitude of control motion to autopilot motion 
against angular frequency. The lower three dashed curves are three 
possible experimental phase-angle curves for the automatic pilot. The 
intersection of the experimental and calculated control-amplitude 
curves establishes the approximate frequency of the airplane with the 
autopilot in operation. If, as in the case of the intermediate experi- 
mental phase-angle curve, the intersection of the experimental and 
calculated phase-angle curves is at the same frequency as the inter- 
section of the control-amplitude-ratio curves, the airplane may be 
neutrally stable and may be expected to osci lla te continuously at this 
frequency. It is, however, more usual that the intersection of the 
experimental and calculated phase-angle curves will not be at the same 
frequency as the intersection of the control-amplitude-ratio curves. 

If the phase-angle curves, as in one case shown, intersect at a higher 
frequency than the control— amplitude— ratio curves, the aircraft will be 
stable. If, as in the remaining case, the intersection of the phase- 
angle curves is at a lower frequency than the control-amplitude— ratio 
curves, the aircraft will be unstable. Because of the nonlinear 
characteristics of the control system, it is generally necessary to 
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make the experiments for different amplitudes. With a dead spot 
(insensitivity to small deviations) there will probably be some ampli- 
tude below which the system will be unstable. 

Calculated frequency— response curves for the aircraft .— In the 
application of the method to the case of aileron control of an aircraft 
or a missile independently stabilized about all three axes, the equation 
of motion for determining the control movement is 

(|t d2 + % i ’)*' c h 8 


The calculated steady— state solution of the aircraft in response to a 

sinusoidal forcing function of unit amplitude 8 = sin cot is 

0 = 0 sin (cut +0). (See reference 1.) The values of 

and 0 are obtained over the desired range of angular frequencies ay 

by the substitution of ia> for D in the equation 
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D 2 + C Zp 
°*8 


D 


This substitution is equivalent to specifying an undamped sinusoidal 
motion and results in the expression A + IB from which can be 

obtained — - \/a 2 + B 2 and 0 = tan -1 The angle 0 may denote 

either a phase lag or lead, depending upon its quadrant. If 0 is in 
the third or fourth quadrant, the control lags behind the displacement, 
but if 0 is in the first or second quadrant, the control leads the 
motion of the airplane. The ratio of the amplitudes 8 euad 0 
is — And may be defined as the control-amplitude ratio K of the 
rmax 

control system, that is, the ratio of maximum control deflection to 
maximum displacement in bank„ A plot of _K and 0 against co shows 
the combination of controlr-amplitude ratio and phase lag or lead neces- 
sary to maintain fixed-amplitude oscillations at any given frequency. 
These res ult s are the calculated frequency— response curves due to a 
sinusoidal motion of the aircraft. 


Determination of equivalent sine wave for the automatic— pilot 
response.— The experimental frequency— response curves are obtained by 
ngri-mR-Hng the automatic pilot sinusoidally at various amplitudes 
through the desired range of frequencies. The control is assumed to 
oscillate at the same frequency as the automatic pilot but because of 
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the physical characteristics of the autopilot the control motion may 
differ widely from a true sine wave and may show arbitrary phase, 
amplitude, or wave— form relations (fig. 2). It is necessary, therefore, 
to determine an equivalent sine— wave response for any arbitrary control 
motion. 

In order to determine the equivalent sine wave for an arbitrary 
control motion, the following relations are assumed: 

(1) The work done per cycle by a control following the equivalent 
sine wave on an aircraft having a harmonic displacement sin cot must 
he the same as that done by the actual control variation. 

(2) The angular Impulse of the equivalent sine wave acting on the 
airplane over a half cycle must equal the change in angular momentum of 
the airplane caused by the actual control motion during the same 
interval. 

The work done by a nonharmonic force &(t) of frequency cu upon 
a harmonic motion sin cot is proportional to 



cos cot dt 


where the period of the oscillation is T = — , and is the coeffi- 

cient of the component cos cut. (See reference 2.) This component of 
the control motion that is out of phase with the aircraft motion is the 
only harmonic of the Fourier series representing the forcing func- 
tion 8(t) which contributes to the work done on the aircraft. The 
flngniflT- impulse is obtained by integrating the curve of control deflec- 
tion against time over a half cycle. This component of the control 
motion in phase with the sinusoidal motion of the aircraft, obtained 
from the second relation, is 


A i ' %lT 6(t) at 


where A^_ is the coefficient of the component sin cut. The condition 
of zero net impulse over a series of cycles may be met by adjusting the 


reference axis for &(t) 


so that 



The control motion 
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may then he expressed as the sum. of the In-phase out-of-phase 
components: 


A]_ sin cot + cos cot 


or 


\f A l 2 + B^ 2 sin (cot + 0) 


The control-amplitude ratio K is equal to wAj_ 2 + Bi 2 and the phase 

—1 ®1 v 

lag of the system 0 is tan — „ The control-amplitude ratio and the 

phase lag or lead are determined ^ram records taken of the osci llat ions 
and plotted against co. In contrast to the calculated frequency- 
response curves which involve the aerodynamic and mass characteristics 
of the aircraft, these experimentally determined curves will he 
functions of the dead spots and various types of lag found in the 
control system. In general, the behavior of the automatic pilot will 
he nonlinear and hence a family of curves showing different phase and 
amplitude relations for different amplitudes of disturbance will he 
obtained. 

Comparison of the calculated frequency-response curves of the 
aircraft and the experimental frequency— response curves of the automatic 
pilot .— The two sets of frequency— response curves show, on the one hand, 

the values of K and 0 necessary for hunting at a given frequency 
and, on the other hand, the actual values of K and 0 obtained 
experimentally at this frequency. In order to determine from these 
curves whether the aircraft will hunt in flight, the following condi- 
tions must be satisfied: 

(1) At a given frequency and amplitude the experimental values 
of K «nd g must agree with the calculated values. 

(2) The motion must be stable for amplitudes larger than the one at 
which the airplane will hunt (as determined from the first condition) . 

The first condition indicates that the control-amplitude ratio and 
phase lag or lead obtained as a result of all types of lag in the 
control system must agree with the combination of K and 0 necessary 
far hunting to exist. The second condition is essential to prevent 
instability if the aircraft is. displaced to amplitudes larger than the 
one at which it will hunt. The aircraft is stable at these larger 
amplitudes if, at the frequency far which the calculated and experi- 
mental control-amplitude ratios are equal, the calculated value of the 
phase leg required for hunting is greater than the experimental value. 
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In other words, the calculated value of 9 is a critical value of the 
lag necessary to cause the aircraft to hunt. If experimental values 
of 9 are less than the critical value, the aircraft motion is damped, 
whereas instability occurs if the experimental value of 9 exceeds the 
calculated critical value (f ig. l) . 

Illustrative case.— The equation of motion in hank of a small 
experimental aircraft tested in the Langley 7— by 10-foot tunnel was 
estimated as 


(0.000245D 2 + O.OO2V5D)0 = -0.26VpS 

Solving the equation for S/0 and substituting io) for D give the 
expression 

5 = 0.000926a) 2 - 0.009260)! 

The resul tant values of K and 0, that is, the calculated frequency- 
response curves for the aircraft, are shown as solid lines in figure 3» 
These curves show that for small values of K, the frequency of the 
steady oscillation is low and the motion will not be sustained unless 
the phase lag is large. As K increases, the frequency of the steady 
osc illation increases but the phase lag required decreases. It is 
important to note that an automatic pilot with a constant time lag 2 
would be unstable at high angular frequencies since the relation between 
angni nr * frequency, phase lag, and time lag is 9 (radians) = 0)2. 

The experimental frequency— response curves were obtained by 
oscillating the automatic pilot at. amplitudes of 10° and 20° through the 
range of desired angular frequencies. The phase lag and control- 
amplitude ratio for the two amplitudes were dete rmine d from records 
similar to figure 2 and are plotted as a function of a) in figure 3 
for two values of control-amplitude ratio K. For this particular 
automatic pilot, the control-amplitude ratio -was independent of ampli- 
tude whereas the phase lag varied with amplitude 0 The results in 
figure 3 indicate that, for each control-amplitude ratio, the experi- 
mental values of 9 are greater than the calculated phase lag, and 
hence the aircraft would be unstable. Unpublished results from wind- 
tunnel tests indicated that the motion was unstable, as predicted from 
the curves of figure 3« 

In an effort to make the aircraft stable, the parameters of the 
automatic pilot were modified and additional wind-tunnel tests were 
performed. The conditions selected for these wind-tunnel tests were, 
however, different from those conditions for which the experimental 
frequency— response curves were obtained and hence no direct prediction 
of the aircraft stability could be made. The results of the wind- 
tunnel tests with the modified automatic pilot indicated a steady 
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oscillation, and records taken of the tests shoved that the values 
of K and 9 agreed very closely with the combination of K and 9 
determined from, the calculated frequency— response curves. The test 
points of figure 4 show the combination of K and 9 far the cases in 
which steady oscillations occurred in the roll tests. 


CALCULATED FREQUENCY RESPONSE OF THE AIRCRAFT 
FOR DAMPED OSCILLATIONS 


The previous calculations of the frequency— response curves were 
based on the assumption that the sinusoidal motion of the aircraft is 
neutrally damped. It is often desirable, however, to determine the 
performance of an automatic control device required to cause the motion 
of the aircraft to damp at a sufficiently rapid rate. Although no 
satisfactory analysis of this problem has been given, a qualitative 
indication of the rate of damping to be expected in a given case may be 
obtained by comparing the measured phase and amplitude, of the control 
to the phase and amplitude calculated to be required to enforce a given 
rate of damping. 


Strictly spe aking , the assumed exponential damping of the motion 
would require an exponential decrease in the response of the autopilot 
at decreasing amplitudes. In general, such a linear response cannot 
be expected and hence the method will require careful Judgment in its 
application. 


The equation of damped motion in bank for the illustrative case 
may be written by adding a real part p to the I m a g inar y root ioi 

[0.000245(-+l + iu>) 2 + 0.00245(-P + ln>)]0 = -£.26455 


where p is given a value as determined by the desired rate of damping. 

0 69^ 

The motion damps to one-half its amplitude in T^ 2 = seconds. 

Solving the equation for 5/0 gives the expression 


§ = 0.000926m 2 - 0.000926p 2 + 0.00926p + i(— 0.00926u> + 0. 00l852a>p) 

0 


The frequency— response curves shown in figure 5 were calculated for 
values of p varying from O.I75 to 8.31. The control-amplitude-ratio 
curves are only plotted for values of p equal to 0 and 5 since they 
are the limiting curves for the Values of p investigated. Figure 5 
indicates that. the control- amplitude ratio is almost independent of p 
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whereas the phase lag decreases as 4 increases. For a control system 
with a given contr ol-ainplitude ratio, therefore, the damping of the 
oscillation increases as the phase lag is reduced. If the oscillation 
is to damp one-half in less than 1 / 7.22 second, the control motion must 
lead the aircraft motion. 

In order to predict quantitatively the stability of the motion of 
an aircraft which damps exponentially, the experimental frequency- 
response curves would have to he obtained for the condition where the 
forced oscillation of the automatic pilot also damps exponentially. 


CONCLUDING REMARKS 


A method for predicting the stability in roll of automatically 
controlled aircraft by a comparison of calculated frequency— response 
curves for the aircraft and experimentally determined frequency- 
response curves for the automatic pilot is presented. The method is 
expected to be useful as a means of establishing the specifications of 
the perfo rmanc e required of the automatic control device for pilotless 
aircraft designed as missiles 0 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., June 13, 19k6 
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Figure 1.- Frequency-response curves for stable, unstable, and neutrally damped motion of the aircraft. 
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Figure 2 0 Typical response of control to sinusoidal oscillation of automatic pilot. 
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Figure 5»— Calculated frequency— response curves of the aircraft for damped oscillations 



